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Abstract 

The primary emphasis of this work on kinetics is to illustrate the a posteri- 
ori approach to applications, where focus on data leads to novel outcomes, rather 
than the a priori tendencies of applied analysis which imposes constructs on the 
nature of the observable. The secondary intention is the development of appropri- 
ate methods consonant with experimental definitions. Chemical kinetic equations 
were largely developed with the assumption of rate constant invariance and in par- 
ticular rate constant determination usually required knowledge of the initial con- 
centrations. These methods could not determine the instantaneous rate constant. 
Previous work based on precise simulation data [ J. Math . Chem 43 (2008) 976- 
1023] of a bidirectional chemical reaction system in equilibrium concluded that the 
rate constants is a function of species concentration through the defined and deter- 
mined reactivity coefficients for at least elementary reactions. Inhomogeneities in 
the reaction medium might also lead to cross-coupling of forces and fluxes, leading 
to concentration dependencies. By focusing on gradients, it is possible to deter- 
mine both the average and instantaneous rate constants that can monitor changes 
in the rate constant with concentration changes as suggested by this theory. Here, 
methods are developed and discussed utilizing nonlinear analysis which does not 
require exact knowledge of initial concentrations. These methods are compared 
with those derived from standard methodology for known chemical reactions stud- 
ied by eminent kineticists and in one case with a reaction whose initial reactant 
concentration was not well determined. These gradient methods are shown to be 
consistent with the ones from standard methods and could readily serve as alterna- 
tives for studies where there are limits or unknowns in the initial conditions, such 
as in the burgeoning fields of astrophysics and astrochemistry, forensics, archeol- 
ogy and biology . All four reactions studied exhibited semi sinusoidal-like change 
with reactant concentration change which standard methods cannot detect, which 
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seems to constitute the observation of a new effect tfiat is not predicted by current 
formulations, where the possibility that the observations are due to artifacts from 
instrumental errors or the optimization method is reasoned as unlikely since the ex- 
periments were conducted by different groups at very different times with different 
classes of reactions. Two broad reasons are given for this observation, and exper- 
iments are suggested that can discriminate between these two effects. Although 
first and second order reactions were investigated here using data from prominent 
experimentalists, the method applies to arbitrary fractional orders by polynomial 
expansion of the rate decay curves where closed form integrated expressions do 
not exist. Integral methods for the above will be investigated next. 

Keywords: numerical nonlinear analysis; orthogonal polynomial expansion; chemical 
reaction rate law; a priori; a posteriori 

1 INTRODUCTION AND METHODS 

Current trends in mathematical appUcations almost always indicate the creation of 
mathematical structures that are then supposed and considered to mirror physical real- 
ity and experimental outcomes. Less common are applications that analyze experimen- 
tal data using as closely as possible the operational definition of variables to elucidate 
the validity or otherwise of theories. The main thrust of this sequel is to illustrate 
research that puts priority on the experimental data as a means to constructing or sug- 
gesting theoretical and mathematical structures. The data from highly empirical field of 
chemical kinetics is used within the scope of the definition of the measured variables to 
discover/uncover new effects; it is suggested here based on the outcome of the analysis 
that mathematical analysis of the data in other fields within the operational definition of 
the empirical variables can elucidate new phenomena and suggest how appropriate and 
consistent theory can be constructed a posteriori, rather than the a priori tendencies 
of many applied mathematicians. One aspect of this culture is the cult of prediction 
and of predictability in the natural sciences where resources are expended in perform- 
ing experiments to verify and substantiate theories. The methods developed here are 
of secondary importance compared to the a posteriori analysis of the data and its out- 
comes; these methods refer to variables which come from the experimental definition. 
For an elementary reaction 

Ai + U2A2 . . . f„^A„^ -> Products (1) 

we define the rate constant k as the factor in the equation 

d[Ai] 

= [A] = -m] (2) 

dt 

where 

i=l 

with lR{t) = [R]{t) in general and lA{t) = [A]i{t) in particular, with the notation 
Iai it)lA2 (*)••• 'a„^ {t) = Iqit). The square brackets denote the concentration of the 
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species, and t is the time parameter. For the above, the order O, which need not be 
integer is defined as O = ^"=1 - Clearly, for the above 

We determine k here directly by various methods of computing average and instanta- 
neous gradients for equation (|2]i. In traditional methods, the integrated rate law expres- 
sion is known for only a handful of integer O values of (|2]i which also require initial 
conditions; no such restrictions apply to the current numerical technique. Another class 
of method is through a least squares optimization of the function R{k) for n datapoints 
defined as 

Rik) = j:(^-kiQmy. (4) 



Then, 



implies 



R'{k) = 0^J2 ^ ''hiU) ) hiU) = 



(5) 



Eq.© does not require iterative methods such as Newton-Raphson's (NR) to determine 
the rate constant. A variant of the R{k) optimization above is found in Sec( l2.0.5T ): 
the reason being that we optimize over an intergrated expression rather than directly 
the rate equation ^ such as (fTsT l for the first order rate constant fciand[T9]for the 
second order constant k2- All variants of the above methods will be discussed in 
sequence in what follows. Most kinetic determinations use logarithmic plots with 
known initial concentrations, although there have been attempts at integral methods 
|[l]|2l|3]|4l|5]|6l|2l|8l|9l[l0l[IT][l2land refs. therein] that dispense with the initial 
concentration. (There may be ambiguities in e.g. |1| concerning choice of indepen- 
dent variables that will be discussed elsewhere). These standard methods all assume 
constancy of the rate constant k, and therefore have not inspired methodology that can 
detect the changes to the rate constant that, according to the detailed results of ref . lfT3l 
sheds important information on the activation energy profile changes due to the force 
fields acting on the reacting species. There are conceivably many other reasons for 
variation with time of the rate constants; they include coupling of inhomogenous tem- 
perature field gradients with chemical species fluxes, leading to physical variable in- 
homogeneities in the reaction cell that modifies the rate of reaction with time. This is 
discussed after the data is presented in what follows. 

There have been detailed and specialized reports and treatises of computational 
techniques over the many decades but these have been sparse and far between. Wiberg 
has [ 14 p. 757] described various more advanced series expansion techniques in con- 
junction with least squares analysis to derive kinetic data. His use of numerical inte- 
gration is confined to solving by Runge-Kutta integration a set of coupled equations, 
such as feature in an enzyme-catalysed reaction lfT4l p. 771]. Wiberg in turn draws upon 
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the collective efforts collated by D. F. Detar lfT5l[T6ll . It seems that Detar's collation 
anticipates to some degree many of the developments cited above in this work's bibliog- 
raphy. A first order treatment of a chemical reaction is given in the program LSKINl 
IfTSl p. 126] requires data and time intervals that are conformable to the Roseveare- 
Guggenheim time interval requirement. LSKIN2 1 16 , p.3] solves for rate constant and 
initial concentrations of a second order reaction based on a series expansion of the in- 
tegrated rate law expression. Here, the curvature would introduce "errors" if a linear 
expansion were used. For both these methods, the constancy of the rate constant fc is a 
basic assumption, which is not the case here. 

Nonlinear analysis (NLA), will be attempted here in preliminary form, in order to 
compute both the instantaneous and average rate constants. We analyze 2 first order 
reactions and one second order one using data from prominent kineticists. In addition, 
we select one first order reaction whose initial concentration index is ambiguous utiliz- 
ing the others as a reference to gauge the likelihood of our result based on NLA; if our 
analysis concurs with the 3 reactions from the literature, then one might be confident 
that the NLA analysis of the ambiguous reaction is reasonably accurate. Important 
experiments in science are conducted under uncontrolled conditions, such as in astro- 
chemical reaction rate determinations and photochemical emission spectra in the Mars 
and Titan atmospheres over the several decades |17, 18, 19|. A similar situation obtains 
in forensic science and archeology and in biological physiological rate determinations. 
The basic methods presented here caters for both controlled and uncontrolled initial 
conditions. 

The 3 first order reactions (i)-(iii) and second order reaction (iv) studied are item- 
ized below: 

(i) the tert butyl chloride hydrolysis reaction in ethanol solvent (80%v/v) at 25°C 
derived from the Year III teaching laboratory of this University (UM) where the 
initial concentration, although determined, is ambiguous. Because of time con- 
straints, the inaccurate Atx3 = 2050/iS cm^^ for (i) was determined by heating 
the reaction vessel at the end of the monitoring to 60°C until there was no ap- 
parent change in the conductivity when equilibrated back at 25°C. Reaction (i) 
involved 0.3mL of the reactant which was dissolved in 50mL of ethanol initially. 
The reaction was conducted at 25°C and monitored over time (minutes) by mea- 
suring conductivity (/iS cm^^) due to the release of H+ and CI ions as shown 
below in (|6]l, 

C4H9C1 + H2O C4H9OH + H+ + cr. (6) 

(ii) the methanolysis of ionized phenyl salicylate derived from the literature |f20), Ta- 
ble 7.1,p.381] with presumably accurate values of both the initial concentration 
and for all data sets of the kinetic run. Reaction (ii) may be written 

PS" +CH3OH ^ MS" +PhOH (7) 
where for the rate law is pseudo first-order expressed as 

rate = fc6[PS]" = /ccpHsOHlpS"]. 
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The methanol concentration is in excess and is effectively constant for the reac- 
tion runs l(20i p.407]. The data for this reaction is given in detail in [20, Table 
7.1], conducted at 30°C where several ionic species are present in the reaction 
solution from KOH, KCl, and H2O electrolytes. 

(iii) the primarily SNl substitution reaction 1251 Table IX,p.2071] of tertiary butyl 
bromide (Bu*Br) with dilute ethyl alcoholic sodium ethoxide in ethanol solvent 
where there concurrently occurs an approximately 20% contribution of an El 
elimination reaction.Reaction (iii) may be written 

SNl 

CHgCBr > Products (8) 

where the solvent was EtOH with initial sodium ethoxide concentration [NaOEt] = 
0.02386N at 25°C. The products consisted of approximately 81% substituted ter- 
tiary butyl ethoxide and 19% olefinic molecules due to El elimination. Hence the 
rate constant here refers to a composite reaction (details in ||23l 2064] and ||25l 
p.2070]. 

(iv) the second order E2 elimination reaction ll23l p. 2059-2060 and Table VII,p.2064] 
with reactants isopropyl bromide (Pr*Br) and sodium ethoxide (NaOEt). Reaction 
(iv) involving isopropyl bromide Pr'Br = (CH3)CHBr(CH3) may be written 

E2 

(CH3)CHBr(CH3) + OEt > CH2CHCH3 + Br + HOEt (9) 

where the isopropyl bromide reacts with the OEt^ ion in EtOH solvent at 25 °C 
to yield 80.3% of the olefinic product with some SN2 substitution with the (OEt) 
functional group fTF Table III,p.2061] according to the kineticists. Further de- 
tails and data appear [23, Table VII,p.2064]. It should be mentioned that the 
E2 reaction was inferred to be second order from prior experimental considera- 
tions since the [NaOEt] concentration reduction from the data exactly coincides 
with the reduction of Pr*Br and was not independently determined. This perhaps 
somewhat experimentally questionable technique may well be the reason for the 
instantaneous rate constant as computed here to be somewhat non-smooth, as will 
be discussed later (see Fig. (12111 for the graph). 

"Units" in the figures and text pertain to the appropriate reaction variable dimen- 
sion, for instance either the absorbance for (i) or the conductivity (/iS cm~^) for (ii) 
below. Either because of evaporation or the temperatures not equilibrating after heat- 
ing, the measured Aoo, it would be inferred that for (i) the measured value is larger 
than the actual one determined from the analysis. Reaction (ii) is very rapid compared 
to (i) and the experimental data plots show high nonlinearity. We denote by A the 
measurement parameter which is the conductivity /iS cm^^ or absorbance A [i20i eqn 
7.24-7.26] for reactions (i) and (ii) respectively; A also refers to the concentration [X] 
of species X for reactions (iii) and (iv). The more accurately determined Aqo = ^00 
for (ii) [20,, Table 7.1,p.381] was at approximately 0.897. Analysis of (ii) give values 
of Aoo = Xoo very close to the experimental ones that suggests that our determina- 
tion for reaction (i) Aoo is correct. The experimental data and number of readings for 
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the determination of rate constants is always related to the method used and the order 
of accuracy required in the study; for Khan |20, Table 7.1,first A column], (reaction 
(ii)) , 14 normal readings over 360 seconds (s) sufficed for Khan's purposes, whereas 
for the practical class (reaction (i)), 36 readings over 55 minutes (mins) were taken. 
The meager 14 readings of (ii) covered a major portion of the nonlinear region of the 
reaction, whereas for (i) the many readings were confined to the near-linear regime. 
Linear proportionality is assumed between A and the extent of reaction x, where the 
first order law (c being the instantaneous concentration, k the general rate constant and 
a the initial concentration) is ^ = —kc = —k{a — x); with Aoo — eta, At — ax and 
A(0) = Ao = axo, integration yields for assumed constant k 

Eqn.dTOb determines k if Ao and Aoo are known. 

The analysis for the latter reactions (ii)-(iv) would provide a reference and indica- 
tion of the predicted value in (i) for the initial concentration, apart from checking for 
overall consistency of the methodology in general situations especially when there is 
doubt concerning the value of the initial concentration. 

The methods presented here applies to any order provided the expressions can be 
expanded as an n-order polynomial of the concentration variable against the time in- 
dependent variable. To get smooth curves that are stable one had to modify and use a 
proper curve-interpolation technique that is stable which does not form sudden kinks 
or points of inflexion and this follows next. 



1.1 Orthogonal polynomial stabilization 

It was discovered that the usual least squares polynomial method using Gaussian elim- 
ination f2T\ Sec.6.2.4,p.318 ]to derive the coefficients of the polynomial was highly 
unstable for npoly > 4, which is a known condition f2T, p.318,Sec 6.2 A]. For higher 
orders, there is in addition the tendency to form kinks and loops in an interpolated 
curve for values between two known intervals. Other methods described in special- 
ized treatises 122] Ch.5, Sec. 5. 7-5. 13] , even if robust and stable,such as the Chebyshev 
approximation required values of the proposed experimental curve at predetermined 
definite points in time, which is outside the control of one using predetermined data 
and so for this work, the least square approximation was stabilized by orthogonal poly- 
nomials |21 , Sec. 6. 3] modified for determination of differentials. It is hoped that the 
method can also be extended to integrals in future investigations. The usual method 
defines the n*'* order polynomial p„ {t) which is then expressed as a sum of square 
terms over the domain of measurement to yield Q in (fTTI) . 

Pn{t) = EjU^J"^^ 

(11) 

Q{f,Pn) = Eill [/i -Pn(^i)]^ 

The Q function is minimized over the polynomial coefficient space. In the orthog- 
onal method adopted here, we express our polynomial expression Pm{t) linearly in 
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coefficients aj of ipj functions that are orthogonal with respect to an inner product def- 
inition. For arbitrary functions /, g, the inner product (/, g) is defined below, together 
with properties of the (fj orthogonal polynomials: 



= Ef=i/(^fe)-.9(ifc) 
Lpj) =0 (i^ j) and {^p^,^p^) ^ 0. 

(pi{t) = {t - bi)ip^^i{t) - Ci(pi^2{t) (i > 1), 

(po{t) — 1, and ipj = for j < 1, 

6,; = {t(pi^i,(pi_i)/{ipi^i,ipi^i) {i>l),bi^O{i <1), 

Ci = (ti^i_i,i^i_2)/(</5i-2,</5i-2) (j > 2),andci = (i < 2). 

We define the m*^ order polynomial and associated aj coefficients as follows: 

Pm{t) = Ejlo^j'^jW 



(12) 



(13) 



(14) 



The recursive definitions for the first and second derivatives are given respectively 

as: 

^[{t) = (p^_i(i)(i-6,) + (^.-i(t)"Q(^U(<)(i > 1) 

(15) 

Here the codes were developed in C/C++ which provides for recursive functions which 
we exploited for the evaluation of all the terms. The experimental data were fitted to 
an TO*'* order expression Am {t) defined below 

m m 

The coefficients hi are all computed recursively, and the derivatives determined 
from ( fTSI ) or from ( fT4l ) and ( fT5] l. Once or Oj are determined, then the gradient to the 
curve A„j (t) is computed as 

m 

The ^Q(t) function of (O is expanded similarly as for Xm{t) for order ra. The or- 
thogonal polynomial method is stable and the mean square error decreases with higher 
polynomial order in general monotonically (where n is used to denote the integer or- 
der), but the differentials are not so stable, because of the contribution of higher order 
coefficients in the differential expression as will be shown. From the form of the of 
the equation that will be developed, the rate constant is determined as the gradient of 
a straight-line graph in the appropriate segment of the graph. However, the curvature 
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of the plot will increase with increasing n, giving a poorer value of k, whereas higher 
values of n would better fit the Xvst curve. Hence inspection of the plots is neces- 
sary to decide on the appropriate n value, where we choose the lowest n value for the 
most linear graph of the expression under consideration that also provides a good X{t) 
fit over a suitable time range over which the k rate constants apply. The orthogonal 
polynomial stabilization method provides good A fits with increasing n, but not gradi- 
ents, so that the onset of sudden changes to the gradient which on physical grounds is 
unreasonable can be used as an indication as which is the best curve to select. There is 
in practice little ambiguity in selecting the appropriate polynomials, as will be demon- 
strated. Reactions (i) and (ii) both gauge initial concentrations in terms of the Aoo ( 
Aoo) or final reading of a physical factor proportional to concentration and the structure 
of the analysis is the same and will therefore be discussed simultaneously, followed by 
reactions (iii) to (iv) , where concentrations are measured directly during the course 
of the reaction, which will be discussed together because the form of the boundary 
conditions and data are of the same class. 

2 ANALYSIS OF REACTIONS (I) and (II) 




10 20 30 40 50 60 
Time/mins 



Figure 1 : Plot of (i) using orthogonal polynomials for various orders n. The the least 
squares deviation goes down dramatically with increasing n, which was found not to 
be the case with the normal non-orthogonal polynomial method. 

Figure ([T]i are plots for the different polynomial orders n for reaction (i). It will 
be noticed that higher n values in general leads to better fits visually; the normal least 
squares method leads to severe kinks and loop formation for ^ n > A which is not 
evident here. The reaction (ii) data covers a far greater domain with respect to half- 
lifetimes with only about 14 points (which is a poor dataset with respect to our methods 
but which still gives quantitatively accurate values); because of the relatively more 
rapid curvature changes, we would expect very different gradient behavior as compared 
to (i) with its stronger linearity. 
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The corresponding plots for reaction (ii) are in Fig. (|2|l. 
1 

0.9 

0.8 

c.0.7 
< 

0.6 
0.5 
0.4- 

100 200 300 400 

Time /s 

Figure 2; Experimental points omitting point at Aoo for reaction (ii) at time ~ 2135s. 
The curve is rather non-linear. 

In view of the nonlinearity, we chose a limited regime to curve fit for polynomial 
order n — 3, 4, 5 in Fig.© and the gradient was computed for the n — 5 polynomial 
to determine the rate constants as it was the only order that gave a smooth curve for 
the first 12 consecutive points in the range; the other orders also gave consistent and 
almost equal gradients except at the extreme end points of the range plotted as depicted 
for example in Figs. (l6l[8lfT2li. 

1 

0.9 

0.8 

^0.7 
< 

0.6 
0.5 
0.4 

50 100 150 200 

Time /s 

Figure 3: Experimental points curve fitted with polynomials of order n = 3, 4, 5. The 
fit for this range is excellent, despite the nonlinear nature of the curve 

UnUke reaction (ii), the Aoo for reaction (i) was ambiguous. The plot of (fTOl l was 
made for the same experimental values with different Aoo's, both higher and lower than 
the supposed experimental value for this reaction. The plots in Fig.Qshows increasing 
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Figure 4: Integrated equationlfTOb plot with Aqo from experiment (0) and from two 
different arbitrary values (1,2) for Aoo, which yields two different values for the rate 
constant due to gradient change. 

ka for decreasing Aoo; the choice Aoo — 1050 leads to a value of ka close to the NLA 
values for the different methods discussed which does not require Aoo, but is also able 
to determine this value by extrapolation. The rate constant from NLA is higher than 
that determined from experiment, implying a lower Aoo value which is consonant with 
evaporation of solvent and/or the non-equilibration of temperature prior to measure- 
ment to determine Aoo. Hence elementary NLA allows one to deduce the accuracy of 
the actual experimental methodology in this example. Except for one section, we shall 
apply NLA based on constant k assumption. We also quote some values of Khan's re- 
sults II20I Table 7.1] in Table ([T]i, where some comment is required. The A absorbance 
is monotonically increasing and at higher time (t) values (see [20, Table 7.1]) the ex- 
perimental A value exceeds the Aoo that is determined by the process of minimizing 
df. Hence the minimization of ^ df with respect to Aoo is taken as a protocol for 
determining the best k value even if it contradicts experimental observation. Further, 
this protocol is highly sensitive to A; a change of 10""^ leads to an approximately ten- 
fold change in k. On the other hand, if ^00 determined from experiment as 0.897 is 
accepted, then then computed rate constant for this value is fc = 2.69 x 10~^s^^ im- 
plying that the uncertainty in k is of the order of ±14 x 10"'^. Hence we can conclude 
that the Khan method is a protocol that accepts as correct the k value that is determined 
by the minimization of A^o for a certain Aoo range (w 0.8980 — .8805), which again 
refers to an unspecified protocol as to the choice of the range. 

The general 1st and 2nd order equations We state the standard integrated forms 
below as a reference that requires specification of initial concentrations in order to 
contrast them to the methods developed here. 
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513.5 


109.4 


8.563 


63.26 


212.7 


227.4 




.8805 


.881 


.882 


.883 


.885 


.887 




19.7 ± .6 


18.1 ± .3 


16.5 ± .1 


15.5 ± .2 


14.2 ± .4 


13.3 ± .5 



Table 1: Some results from reaction (ii) ll20l p.381. Table 7. 11. The first row refers to the 
square difference summed, where the lowest value would in principle refer to the most 
accurate value (third entry from left). The second row refers to the Aoo absorbance and 
the last to the corresponding rate constant with the most accurate believed at the stated 
units to be at 16.5 ± .1. 



The first order rate constant ki is determined from 



fci = ilogJ6(6-x)} 



(18) 



and the second order rate constant k2 is determined from 



t{a - b) 



Ids 



b{a — x) 



(19) 



where these expressions are given in 11231 p. 2063] and utilized to compute rate constants 
where a and b are the initial concentration terms at i = and x is the extent of reaction. 



2.0.1 Method 1 

This method is a variant of the direct method of eq.©. For constant k, the rate equation 

^ — —kc — —k{a — x) reduces to 

^ = -kX{t) + Xoo.k (20) 

Hence a plot of vs X{t) would be linear. We find this to be the case for polynomial 
order npoly < 3 as in Fig.© below for all data values; higher polynomial orders can 
be used in selected data points of the curve below, especially in the central region. 
Thus criteria must be set up to determine the appropriate regime of data points for a 
particular polynomial order in NLA. For n=2, ka = 3.34 ± .03 x 10^'^min^^ and 
Aoo — 1134 ± 10 units. The plots for reaction (ii) is a little more involved; it is a much 
more rapid reaction and the number of data-points are relatively sparse for NLA and the 
points cover the entire range of the reaction sequence and is highly non-linear; it was 
found that the gradients were smooth for the first 10 or so points and reasonably linear, 
but that at the boundary of these selected points, there are deflections in the curve; on 
the other hand , the different order polynomial curves (n < 5) are all coincident over 
a significant range of these values; we chose the n — 5 polynomial curve to determine 
the curve over the entire range and the linear least squares fit yields the following data 
kb = 1.64 ± .04 X lO^^g-i and = 0.8787 ± .0008 units. 
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Figure 5: Method 1 graph showing hnearity lower order polynomial fits for reaction 
(i). 




Figure 6: Method 1 apphed to reaction (ii). Only at the peripheral value does the fit fail 
for lower values of n due to the extreme curvature. 
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2.0.2 Method 2 

This method is yet another variant of the direct method of eq.(|2|l. Let a' = Xoo — Xq, 
then In a' — ln(Aoo — X) = kt, then noting this and differentiating yields 




Time /mins 



Figure 7: Method 2 reaction (i) where smooth curves are obtained for at least npoly < 
4. 




Figure 8: Method 2 reaction (ii) where smooth curves are obtained for npoly < 5. The 
n — 5 polynomial is used to calculate the best linear line over this range as the other 
polynomials do not fit for the last value in the series due to the extreme curvature. 
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A typical plot that can extract fc as a linear plot of lii(c?A/di) vs t is given in Fig.© 
for Method 2, reaction (i) and in Fig.® for Method 2, reaction (ii). Linearity is ob- 
served for npoly = 2 and smooth curves without oscillations for at least npoly < 3 
for reaction (i) and an analysis for reaction (ii) uses npoly = 5. The linear least square 
line yields for Method 2 the following: 

ka = 3.35 ± .03 X lO^^min"^ and Aoo = 1130 ± 10 units 
kb = 1.72 ± .02 X lO^^g-i and Aoo = 0.86(53) ± .02 units. 

We note that because of the manifest nonlinearly of the gradients, one cannot de- 
termine the Aoo values to 4-decimal place accuracy as quoted by Khan based on his 
model and assumptions |20, Table 7.1]. 



2.0.3 Other associated non-direct methods 

There are other methods, one of which is a variant of Method 1 and another that utilizes 
a least-squares optimization of the form of the equation for first and second derivatives. 



2.0.4 Method 1 variation 

A variant method similar to the Guggenheim method ^ of elimination is given below 
but where gradients to the conductivity curve is required, and where the average over 
all pairs is required; the equation follows from (ISTT i. 

N N 

(fc) = ^^^7^^r3ITEEln(^'(^«)/^'fe))/(^* -^^O (22) 

^ ' i j>i 

Since we are averaging over instantaneous k values, there would be a noticeable stan- 
dard deviation in the results if the hypothesis of change of rate constant with species 
concentration is correct. Differentiating ( |2TI ) for constant k leads to ( |23] ) expressed in 
two ways 

d^A , fdX\ , , , £X , fdX\ 

. (a)orfc^— — / — (6) (23) 



dt^ \dt J ' dt^ \dt J 

If A(i) = X]r=!o^ a{i)t^^^, then as < ^ 0, the rate constant is given by fc = ^^"j^^^-' 
from (l23b). For the above, n,id, andm denotes as usual the polynomial degree,the 
lower coordinate index and the upper index of consecutive coordinate points respec- 
tively, where the average is over the consecutive points, whereas the fc rate constant 
with subscript "all" below refers to the equation (l22l i . 

The results from this calculation are as follows: 

fca,ai/,fca,zd,z« = 3.32, 3.23 ±.07 X 10-^min-^ n = 2,id= 10,zu = 20 

ka,t^o = 3.082 X lO-^min^^ 
kb,aii,kb,id,iu — 1.7150, 1.676 ± .3 x 10~^, s~^, n = 5,id = l,iu = 10 
kbj-,a = 1.023 X lO-^s"^ 

The asymptotic limit gives a lower value for fcf, than for the other methods for 
reactions (i) and (ii). One possible explanation is that the rate constant changes as a 
function of time, but we note that ( |23] ) was derived assuming constant fc. 
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2.0.5 Optimization of first and second derivative expressions 

Eq.(l23l'b)) suggests another way of computing k for "well-behaved" values of the dif- 
ferentials, meaning regions where k would appear to be a reasonable constant. The (a) 
form suggests an exponential solution. Define ^ = dl and = d2l. Then dl{t) = 
Aexp{—kt) and dl{0) ^ A = h2 from (fTSIl . Furthermore . as t 0, k — (— 2/12//11) 
and a global definition of the rate constant becomes possible based on the total system 
A(i) curve. 

With a slight change of notation, we now define dl and d2l as referring to the 
continuous functions dl{t) ~ Acxp{—kt) and d2l{t) ~ -'kAexp{—kt) and we con- 
sider (dX/dt) and d^X/dt^ to belong to the values (fT6b derived from Is fitting where 
(dX/dt) = A'„, [d^X/dt^) — A'^ which are the experimental values for a curve fit 
of order m. From the experimentally derived gradients and differentials, we can define 
two non-negative functions Ra{k) and Rp{k) as below: 
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R^ik) = - dl{U)) (24) 

where 

Uik)^R'^{k) and Jp{k) = R'^{k) 

and a stationary point (minimum) exists at /Q.(fc) — fp{k) — 0. We solve the equations 
fa , fp for their roots in k using the Newton-Raphson method to compute the roots 
as the rate constants ka and kp for functions fa{k) and //^(fc) respectively. The error 
threshold in the Newton-Raphson method was set at e = 1.0 x 10^^ We provide a 
series of data of the form [n, A, k^, kp, Xa.oo, ^p.oo] where n refers to the polynomial 
degree, A the initial value constant as above, ka and kp are the rate constants for the 
functions /„ and fp (solved when the functions are zero respectively ) and likewise for 
Aq.oo and X/s^oo- The e symbol refers to base 10 (decimal) exponents. The Aoo values 
are averaged over all the 36 data points for reaction (i) and for the 12 datapoints of 
reaction (ii) from the equation 

X^^^Uxit) (25) 
dt k 

for scheme a and /3 for both reaction (i) and (ii). The results are as follows. 
Reaction (i) 

[2, 3.7632 X 10°, 3.2876 X lO^^, 3.2967 x lO^^, 1.1506 x 10^, 1.1477 x 10^], 
[3, 3.6384 X 10°, 2.7537 X lO^^, 2.7849 x lO^^, 1.34756 x 10^, 1.3334 x 10^], 
[4, 3.6384 X 10°, 2.0973 X 10"^ 2.4716 x 10"^ 1.7408 x lO^, 1.4900 x 10^], 
[5,4.0213 X 10°, 9.7622 x 10-^,4.9932 x 10^3,4.4709 x 10^,7.9328 x 10^,], 
[6, 4.5260 X 10°, 4.1270 x 10"^ 8.9257 x lO^^, 1.7101 x 10^, 4.8403 x 10^] . 
We noticed as in the previous cases that the most linear values occur for 1 < rt < 4. In 
this approach, we can use the and fp function similarity of solution for k^ and A;^ to 
determine the appropriate regime for a reasonable solution. Here, we notice a sudden 
departure of similar value between ka and kp (about 0.4 difference ) at n = 4 and so 



15 



we conclude that the probable average "rate constant" is about the range given by the 
values spanning n = 2 and n ~ 3. Interestingly, the Aqo values are approximately sim- 
ilar to the ones for method 1 and 2 for polynomial order 2 and 3 for reaction (i). More 
study with reliable data needs to be done in order to discern and select appropriate 
criteria that can be applied to these non-linear methods. Because of the large number 
of datapoints in the linear range, ka and kt values are very compatible for n — 2,3 
where the ka determination involves double derivatives, which cannot be determined 
with accuracy unless a sufficient number of points is used. 
Reaction (ii) 

The results for this system are 

[5, 7.5045 X 10-3, 1.2855 x lO^^ 1.5497 x lO^^ .94352, .89247] 
for the first 12 datapoints of the published data to time coordinate 155 sees. For polyno- 
mial order 3,4 and the first 1 1 datapoints, where there are no singularities in the curve 
we have 

[3, 7.7275 x 10-3, 1.4469 x IQ-^, 1.6147 x 10-^, .91320, .88335] 
[4, 7.4989 X 10-3, 1.3146 x 10~2, 1.5359 x 10-^, .94208, .89652]. 
Here, ka and ki, differ by ^ .2 x lO^s"^; one possible reason for this discrepancy is 
the insufficient number of datapoints to to accurately determine Even if the num- 
ber of points are large, experimental fluctuations would induce changes in the second 
derivative which would be one reason for discrepancies. Hence experimentalists who 
wish to employ NLA must provide more experimental points, especially at the linear 
region of the X{t) vs t curve. 

2.1 Inverse Calculation 

Rarely are experimental curves compared with the ones that must obtain from the ki- 
netic calculations. Since the kinetic data is the ultimate basis for deciding on values of 
the kinetic parameters, replotting the curves with the calculated parameters to obtain 
the most fitting curve to experiment would serve as one method to determine the best 
method amongst several. For reaction (i) we have the following data: 



Result 


Procedure 


Poly, order 






1 


From expt 




2050. 


1.7752 X 10-3 


2 


Method 1 


2 


1134.3 


3.3397 X 10-3 


3 


Method 2 


2 


1130.23 


3.347 X 10-3 


4 


sec(l2.0.5b i?« 


2 


1150.63 


3.288 X 10-3 



Table 2: Data for the plot of Fig.® for reaction(i) 



For reaction(ii), we have the following: 

Fig.® indicate that the parameters derived from experiment is the worst fit com- 
pared to the methods developed here, verifying that our computations, including the 
Aoo values are a better fit than the one derived from experiment due to flaws in the 
methodology of driving the reaction to completion by heating, leading to evaporation 
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Result 


Procedure 


Poly, order 






1 


From expt 




.8820 


1.65 X 10"^ 


2 


Method 1 


5 


.8787 


1.64 X 10"^ 


3 


Method 2 


5 


.8653 


1.72 X 10"^ 


4 


sec(l2.0.51l Rr 


5 


.89247 


1.5497 X 10-^ 



Table 3: Data for the plot of Fig.dTOll for reaction (ii). 
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Figure 9: The plots according to the ki, and Aiufty values of Table© 
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Figure 10: The plots according to the fcf, and Aiufty values of Table© 



and therefore inaccurate determination. Based on the comparisons between the reac- 
tions (i) and (ii), we predict that reaction (i) if carried out under stringently controlled 
conditions, especially in determining Aoo would have a rate constant approximately 
~ 3.2 X 10^'^mins^^ rather than the experimentally deduced ^ 1.77 x lO^'^mins^^ 
with Aoo ^1130 units rather than 2050 units . For reaction (ii), we note a good fit for 
all the curves, that of the experiment. Khan's results and ours. 



2.2 Evidence of varying kinetic coefficient k for reactions (i) and 
(ii) 

Finally, what of direct methods that do not assume the constancy of k which was the 
case in the above subsections? Under the linearity assumption x ~ a\{t), the rate law 
has the form dc/dt — —k{t){a — x) where k{t) is the instantaneous rate constant and 
this form implies 

If Aoo is known from accurate experiments or from our computed estimates, then k{t) 
is determined; the variation of k{t) could provide crucial information concerning reac- 
tion kinetic mechanism and energetics, from at least one theory recently developed for 
elementary reactions I.13J at equilibrium; and for such similar theories [24| and experi- 
mental developments for very large changes in concentration, it may be anticipated that 
nonlinear methods would be used to accurately determine k{t) that would yield the so- 
called "reactivity coefficients" 1 131 that account for variations in k that would provide 
fundamental information concerning activation and free energy changes. However, 
since these coefficients pertain to the steady-state scenario where a precise relationship 
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exists between the ratio of these coefficients and that of the activity, one might not ex- 
pect to detect these coefficients for relatively minute concentration changes that occur 
in most routine chemical reaction determinations. On the other hand, the very prelim- 
inary results here seem to indicate transient variations belonging possibly to another 
class of phenomena ; it could well be due to periodicity in the reactions where some 
of the "beats" detected -assuming no experimental error in the data- could be related to 
the time interval between measurements, that is, because the number of datapoints is 
restricted, only certain beats are observed in the periodicity. The assumption that the 
spectroscopic detector is not noisy relative to the magnitude of the experimental data 
and that it does not have significant periodic drift relative to a reference absorbance 
leads to the conclusion that some form of chemically induced periodicity might be 
present. It would be very interesting to increase the number of datapoints where the 
time interval between measurements is reduced and to analyze the different types of 
apparent frequencies that might be observed with different time intervals of measure- 
ment. From these observations, perhaps theories could be adduced on the nature of 
these presumed fluctuations. 
Reaction (i) results 

FigurelfTTTi refers to the computations under the assumption of first order linearity of 
concentration and the conductivity. Whilst very preliminary, non-constancy of the rate 
constants are evident, and one can therefore expect that another area of fruitful exper- 
imental and theoretical development can be expected from these results that incorpo- 
rates at least these two effects. 
Reaction (ii) results 
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Figure 1 1 : Variation of k with time or concentration changes based on the experimental 
value Aoo = 2050 units and the computations based on different polynomial degrees 
n = 2, 3, 4 and the computed Aoo values for Method 1 and Method 2 for fixed polyno- 
mial degree n — ?>. 

To verify that the curious results are not due to minute fluctuations of the gradient. 
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we plot the gradient ^ of the curve fits for polynomial orders 1,2,4 and 5. Even for 
low orders, the fit is very good with no oscillations observed between lower and higher 
order polynomials for approximately the first 10 values of the kinetic data in Fig. (fT2l l. 
It was found that the lower order polynomials gave essentially the same results for the 
restricted domain where the gradients coincided with those of higher order 
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Figure 12: Variation of gradient dX = dt = dA/dt with time for different polynomial 
orders n — 1,2,4 and 5 

The gradient drops to at the long time i — > oo limit; on the other hand, the 
factor -jj — -\{t)) ^^^^^ to infinity; so we might expect from these two competing factors 
various sinusoidal-like properties, or even maxima. The surprising result is shown in 
Fig. ( fT3] l. It could be that the form ( l26b is not valid because no instantaneous value 
of the rate constant can be defined. Also, it is not possible at this stage to definitively 
rule out the detector causing a sinusoidal variation due to periodic drift and instability. 
On the other hand, if we can rule out artifacts due to systematic instrument error, then 
one must admit the possibility that a long-time cooperative effect involving coupling 
of the reactant molecules through the solvent matrix over the entire reaction chamber 
may take place . This seems like a big idea that could be investigated further. However, 
if the results are due solely to instrument error, then this method allows us to monitor 
the error fluctuations by taking a suitable weighted average. We note that reaction (ii) 
is relatively complex, involving many ionic species and some intermediate steps or 
reactions 1*20^, p. 414-416]. This fact, coupled with the cell setup where steady state 
temperature gradients might well exist, would lead to coupled processes described by 
irreversible thermodynamics, which could possibly explain such rate constant changes 
relative to the first order parametrization used here. 

Comment: Barring artifacts, Figs. (fTT]fT3T l is consonant with two separate effects: 
(a) a long-time Umit due to changes in concentration that alters the force fields and 
consequently the mean rate constant value (according to the theories in ifTSl l24l ) of 
the reaction as equilibrium is reached, and (f3) possible transient effects due to collec- 
tive modes of the coupling between the reacting molecules and the bulk solution as 
observed in the region between the start of the reaction and the long-time interval. In 
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Figure 13: An apparently periodically varying rate constant that settles to a higher 
value at larger time increments. Aoo was the value taken from Method 1 above for 
polynomial order 5. 

both reactions (i) and (ii), there appears a slight change in the rate constant value at 
time t = and the values at the end of the experimental measured interval which may 
be due to the altered force fields that would change slightly the rate constant according 
to (a). On the other hand, there is a relatively slow and minute sinusoidal-like change 
in the rate constant that may be due to some cooperative effect, if no artifacts are im- 
plied; the interpolation with different polynomials leading to the same gradient seems 
to suggest that some type of collective behavior might be operating during the course 
of the reaction; if this is so then (/?) would be a new type of phenomena that has not 
hitherto been incorporated into chemical kinetics research. 

3 RESULTS FOR REACTION (III) and (IV) 

Two different methods are utilized to determine the reaction rate constants. The direct 
method utilizes determining the gradient k of the d[Ai]/dtvs [Q] curve of eq.(|2]i by 
fitting the best straight line. Initial concentrations are not required, and the error in 
the gradient may be estimated from the mean least squares error of the end-points; 
define the mean square per point A"^ as A'^ = ~ Vif' /"^ where y is the linear 

optimized curve and a datapoint within the range of measurement. Then for the range 
of datapoints \X\ = \Q\ we estimate the error in the rate constant as A/c = \/A^/|X|. 
For what follows below the first order rate constant for reaction (iii) is denoted fcid 
derived from direct computation of the gradient of (|2| whereas kus denotes the rate 
constant as calculated by our least squares method Q. Similarly, for reaction (iv), 
k2d denotes the rate constant derived directly from the gradient of the curve following 
eq.© by fitting the best straight line and k2is is the second order rate constant from the 
least squares minimization of (|5]). A detailed description follows below. 
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3.1 Reaction (iii) first order details 

Fig.(fT4li is a plot of the various orthogonal polynomial order fits. The n — 2 order is 
rather poor but higher orders all coincide with the experimental points. And for n > 2, 
we find that the gradient curves coincide within a certain range where [A] > 0.005M. 
If in fact the order is 1 or unity, a plot of [A]t)s[A] would be entirely linear. Fig. dTsl l 
depict these plots, where some linearity is observed for [A] > 0.005M. 

0.01 8| ■ . ■ ■ . 1 




Figure 14: The polynomial degree n for the orthogonal polynomial fit for the first order 
Bu*Cl reaction (iii) where there is near coincidence for n = 4, 5, 6, 7. 

Fig.dTSIl shows a coincidence of the rate curves for order n — above [A]=0.005M. 
It may be therefore inferred that for at least n > 5, and for [A] > 0.005M, the gradi- 
ent represents the rate constant. The inaccuracy for very low concentrations may be 
explained by reference to Fig.(fl4]i. The experimental curve is parametrized as 

[A] = J2 ^^t' (27) 

i=l 

The time parameter is very large at low concentrations (> 3.5 x lO^s) of Bu*Br; the 
hi values would be small and the uncertainties in the reactant concentration relatively 
high; this explains the large scatter in the gradient values at low concentrations. We 
therefore ignore the first two values of [A] at low concentrations and focus on the 
gradients for points of coincidence of the different polynomial curves in in Fig.lfTsTl. 
The results are shown in Table |4] The linear fit in this specified range yields fci^. The 
absolute root mean squared deviation per datapoint of Fig.(fT4li is listed in the last rhs 
column of Table|4]where the best fit is in the range n = 5 — 7. The average ki^ value 
taking into account the error estimate is (5.4 ± 0.5) x lO^^s^^. which is close to the 
(5.22 ± 0.3) X lO^^s^^ of the experimentalists. For the same regime 

The results from the literature for this first order reaction 1251 Table IX,p.2071] us- 
ing equation (fTSl l has a mean value of 5.22 x 10~^ and for the 12 datapoints determined 
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Figure 15: A plot of the rate of decomposition of Bu*Br(= A)vsA where [A] refers to 
the molar concentration of Bu*Br according to the kinetic data published in ll25l Table 
IX.] for reaction (iii). 



n 


fcid/s"^ 


fcl/s/s-^ 


Est. error in k 


Abs. Dev. in n 


2 


5.038689x10-6 


4.881742x10-6 


4.164406x10-^ 


0.000151 


3 


5.374684x10-6 


5.307609x10-6 


8.129860x10-8 


0.000073 


4 


5.053045x10-6 


5.047342x10-6 


3.541633x10-^ 


0.000034 


5 


5.398782x10-6 


5.181527x10-6 


2.772767 xlO-'^ 


0.000014 


6 


5.393767x10-6 


5.184763x10-6 


2.723767x10-^ 


0.000013 


7 


5.500610x10-6 


5.187774x10-6 


3.028407 xlO-'^ 


0.000014 



Table 4: Results for the first order Bu*Br reaction neglecting last datapoint for calcu- 
lating the mean rate constant k. 



, and this value of k varied with range (5.04 — 5.48) x 10-6 in appropriate dimensions. 
The results as computed according our methods in Table |4] Table |4] shows the absolute 
deviation per point to be quite small, and Fig.lO shows some plots. The gradient of 
the curves are graphed for various n in Fig.(fT5]l. The linear fit in this specified range 
yields kid- The absolute root mean squared deviation per datapoint of Fig. (fT4b is listed 
in the last rhs column of Table |4] where the best fit is in the range n = 5 — 7. The aver- 
age kid value taking into account the error estimate is (5.4 ± 0.5) x 10-6s-^ which is 
close to the (5.22 ± 0.3) x 10-6s-^ of the experimentalists. For the same regime, (|5]l 
is used to compute kus listed in the 3rd column. For n — 5 — 7, kis = 5.18 x 10-6s-^ 
which is exceptionally close to the experimental determination mentioned above based 
on the integrated equation ( fTSl l. requiring initial concentrations. 

Lastly, (O is used to compute the instantaneous rate constant shown in Fig.(fT6]l. 
We note that a maximum is formed before a drop at lower concentrations. Again, the 
concave form with a maximum is evident here as for reactions (i)-(ii). 
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Figure 16: The direct calculation of the change of the rate constant with concentration 

[A] =[Bu*Br] directly from the published experimental data Il25l Table IX.] for reaction 
(iii). 

3.2 Reaction (iv) second order details 

The pioneer experimentalists had decided a priori that the NaOEt iso-PropylBr (Pr'Br) 
was second order and therefore did not independently measure the NaOEt and the Pr'Br 
concentrations as shown in Fig.lfTSI). Such a setting introduces a larger degree of scat- 
tering even if the rate order in known a priori to be the appropriate one; the scattering 
becomes evident in the rate constant curve of Fig.(l2Tli. The orthogonal polynomial fit 
is very good for ?i > 3 in Fig.dTSb. The gradient curve of Fig.(fT9l) shows a coincidence 
of points for polynomial order n > 3 except for the measurement at the lowest con- 
centrations, for the same reasons as given for reaction (iii). Fig. ( l20b is a close-up of 
the rate vs [A] [B] curve where the first 2 points show significant scatter. We ignore the 
first 3 points of lowest concentrations in our calculations for k2d and k2is for different 
polynomial orders n; k2d is the rate constant by linear least squares fit to each of the 
curves of Fig. ( fT9] l of d[A]/dt vs [A][B] for different polynomial orders n and k2is is 
the rate constant calculated according to (|5]l. The results are presented in Table |5] In 
Table |5] the value of k2d is remarkably constant for ?i = 4 to 7, in keeping with the 
curves of Fig. iT% that are coincident for the selected concentration ranges where we 
have k2d = (2.80 ± .07) x lO^^M^^s^^ which is very close also to the computed k2is 
values, where the average value may be written k2is — (2.83 ± .05) x 10~^M^^s^^. 

The experimental results for this second order reaction [23 Table Vll,p.2064] using 
equation (fT9]l has a mean value of 2.88 x lO^^M-^s^i (2.95 x lO^^M^is^i if solvent 
expansion is taken into account according to an extraneous theory), and for the 19 dat- 
apoints determined , k varied within the range (2.76 — 2.96) x 10~®)M~^s~^. Fig. (12111 
graphs the instantaneous rate constant for reaction (iv) for for the concentration ranges 
used for calculating k2d and k2is- Again a semi-sinusoidal shaped curve is observed. 
There is evident scatter in this graph which may be attributed to the fact that [A] and 

[B] are completely correlated and there is no possibihty of random cancellations due to 
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independent measurement of both [A] and [B]. The form of the curve is as for reaction 
(iii) in Fig.(fT6ll and as for reaction (ii) in Fig.(fT3b if we ignore the low concentration 
value as having a large scatter at about t — 150s. Even the ambiguous reaction (i) 
shows a shallow concave shape, as with all the rest. Fig.(fT7]i is a pallet for the 4 diverse 
reactions (i)-(iv) of different orders reported by different sources all depicting the same 
general form; suggestive of some type of overall "chemical inertia" effect. The data 
for reaction (ii) were determined over 50 years apart from (iii) and (iv) by prominent 
persons, especially Ingold and co-workers who were amongst the best kineticists of the 
20th century, apart from elucidating and defining the SNl, SN2, El and E2 reactions in 
organic chemistry. It seems that all these reactions depict a transient and long-ranging 
coupling phenomena hitherto unnoticed due to traditional analysis that uses integrated 
rate law expressions with the presupposition of invariant rate constants, which is also 
built into current statistical mechanics theories. There could be at least two separate 
possibilities: 

(a) independently of the initial concentrations, one observes a rise in the rate con- 
stant before falling when the reactant concentration falls. Hence if we started a ractin 
at the midpoint concentrations close to the peak rate constant value in the reaction runs 
for reactions (ii)-(iv), then the rate constant profile will show a form similar to those 
shown here, except the peak rate constant will shift to lower concentrations, at least 
at a concentration less than the concentration at the commencement of the reaction. 
This possibility suggests some type of "chemical momentum" which is a long-ranging 
coupling phenomena that current statistical mechanical theories are not able to account 
for. Obviously this scenario admits the possibility of hysteresis behavior 

(b) the rate constant is simply a function of reactant and product concentrations as 
outlined for instance in 1 13| and no hysteresis behavior exists. 

Currently, mainstream statistical mechanical theories do not have a quantum or 
classical description of (a) and (b) above. 

Reality would probably be described by broadly either (a) or (b) above with a pos- 
sibility of the other playing a minor role in the contributed effect. 

More careful experiments under stringent conditions need to be performed to : 

a. verify the existence of this effect by using other methods just in case the polynomial 
method has a property that induces a maximum in the gradient at approximately 
the midpoint of the domain range under investigation 

/3. determine which of the two (a) or (b) above is the preponderant effect. 

4 CONCLUSIONS 

The differential methods developed here yield results for the average rate constant that 
is consistent and close in value to the traditional integrated law expression for 3 re- 
actions whose rate constants were determined with precision by reliable and promi- 
nent kineticists, implying that the methods developed here are robust. Based on re- 
producibility of the method in relation to the standard protocol, we further test the 
differential methods for the ambiguous reaction (i) with regard to the inaccurate to- 
tal reactant concentration (reflected in Aoo) and we showed that the rate constant can 
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(c) Reaction (iii) (d) Reaction (iv) 



Figure 17: Collation of change of rate constants with concentration and time for reac- 
tions (i) to (iv). 



n 


fc2d/M-is-i 


Error Afc2d 


A;2;,s/M-is-i 


2 


1.251043x10-*^ 


1.189216x10-^ 


2.180891x10-6 


3 


2.142702x10-6 


8.845014x10-* 


2.654936x10-6 


4 


2.818337x10-6 


6.850636x10-^ 


2.829915x10-6 


5 


2.822939x10"" 


6.600369x10-'' 


2.833440x10-6 


6 


2.876573x10-6 


6.613089x10-* 


2.853618x10-6 


7 


2.770475x10-6 


7.707235x10-* 


2.823186x10-6 



Table 5: Results for the second order Pr'Br- NaOEt reaction neglecting the last 3 data- 
points that are discontinuous for calculating the rate constants k2d and k2is- 
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Figure 18: The experimental decay curve for the NaOEt iso-PropylBr (Pr'Br) reaction 
123. Table Vll,p.2064] for reaction(iv). 
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Figure 19: The rate curve shows that the polynomial order n=2-3 is to low to fit the 
changes of gradient. The other curves of higher order all coincide exact for very large 
time values or low reactant concentrations. The average gradient is calculated by ig- 
noring these values that are off-scale. [[A] = [NaOEt]] and [[B] = [iso-PropylBr]] in 
Molar concentration units. 
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Figure 20: Close-up of the rate curve with reactant product concentration. The last 3 
lowest product concentration points are ignored in the rate averaging to determine the 
rate constant. 




Figure 2 1 : The computation of the instantaneous rate along the regime of coincidence 
of the polynomials for n = 4 — 7 forn = 6. 
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be determined without such data. Since much science refers to systems whose key 
variables are not controllable, as in astrophysical measurements and in forensics and 
archaeological studies, these methods could prove useful for analysis in these areas for 
reactions to arbitrary order. 

All the reactions studied with different orders and mechanisms within the polyno- 
mial optimization method all show some type of "chemical rate constant momentum" 
effect in that there is a gradual acceleration in the rate over time initially followed by 
a sharp decrease as reactants depleted. This observation appears to be novel. Two 
independent factors might contribute to this effect: (A.) the transient coupling of ther- 
modynamical force gradients with flows and species concentrations and molecular ori- 
entation that leads to the observed profile 

(B.) the rate constant is a second order function of the reactant and product con- 
centration according to [ 13 1. That theory was based on evidence from equilibrium MD 
simulation of a chemical reaction. 

Stringently controlled repeated experimentation is required to determine the rela- 
tive contributions of (A.) and (B.) above. If a semi-sinusoidal profile is observed for 
different initial concentrations that covers the range of the concentration profile of a ref- 
erence reaction of exactly the same type which also exhibits a semi-sinusoidal profile, 
then (A.) is the predominant mechanism which exhibits hysteresis behavior, whereas if 
the resulting profile of the experiments leads to a truncated semi-sinusoidal curve of the 
reference profile beginning at the initial concentration of the experimental run, the (B.) 
is the major effect. There is of course the possibility of a combination of effects (A.) 
and (B.) to varying degrees. Both effects have not been anticipated in current kinetic 
theories, which implies that these effects in themselves constitutes one area for further 
investigation. 

Whilst the main purpose of this work is not to provide physical explanations, we 
suggest that the resulting curves in Figs.(fT7]) can be explained by introducing X con- 
tributing factors or processes in real time. The first (fl) involves the fact that the re- 
actants are initially separated, and the molecules must diffuse homogeneously before 
they can begin to interact. The second process (f2) involves reactant interaction with 
the solvent matrix, which would impede the reactive interactions and also conceivably 
raise the activation energy relative to unbounded reactants and lastly (f3) describes the 
product solvent matrix interaction, which is probably is not too significant for the typ- 
ical reactions studied here. At the initial stage of the reaction, there is minimal solvent 
reactant interaction which would result in the caging of the reaction active sites and 
so the reaction is diffusion limiting; within a certain time scale, the mutual diffusion 
of reactants would allow for more reactant-reactant interactions, leading to an apparent 
increase in the rate constant; with the progress in time, however, the caging of reactants 
due to (f2) would increase the effective activation energy and hence lower the value of 
the rate constant which explains the precipitous drop at large time values; process (f3) 
might prevent in some cases the back reaction due to the breakup of the product to 
reactant molecules, and it may moderate (f2) by liming the number of active solvent 
interaction with the reactant molecules. One might expect (fl) to case the rise in the 
rate constant, and (f2) the lowering, leading to the concave maximum observed in all 
the reactions due to these competing processes. Conventional kinetics, modeled af- 
ter ideal situations of homogeneity, is not able to account for these fine second order 
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details in the change in the reaction rate. Even in the homogeneous case , such as a 
reactive system in thermodynamical equihbrium, it was found that the rate constant 
is a well-defined function of the reactant and product concentrations, but this contri- 
bution to the changes found here is probably of second order for the typical reactions 
mentioned here. 

The results presented here provides alternative developments based on NLA that 
is able to probe into the finer details of kinetic phenomena than what the standard 
representations allow for, especially in the the areas of changes of the rate constant 
with the reaction environment as well as determine average rate constant values without 
Aoobeing known. Even with the assumption of invariance of k, one can always choose 
the best type of polynomial order that is consistent with the assumption, and it appears 
that the initial concentration as well as the rate constant seems be be predicted as global 
properties based on the polynomial expansion.lt should be noted that the examples 
chosen here were first order ones; the methods are general and they pertain to any 
form of rate law where the gradients and forms can be curve-fitted and the form of 
the equations optimized as in section ( I2.0.5I I. One other research area that may be 
investigated is the possibility of reactions of fractional order; elementary reactions are 
by nature of integer order; is there a method that can reduce them to fractional order if 
the rate constant is indeed in part a weak function of the reactant concentrations? 
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